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Abstract 

The Cornell potential with the best fitted parameters a* and k are modified by adding terms 
derived from non-perturbative QCD, which are characterized by a series of non-vanishing 
vacuum condensates of quarks and gluons. In terms of this potential, we study the system of 
heavy quarkonia. The results show that the correction caused by the additional terms reduces 
the deviation between the data and the values calculated with the pure Cornell potential and 
improves the splittings of energy levels. The achievements indicate that the non-perturbative 
effects induced by vacuum condensates play an important role for the correction to 1/q^, 
which in general was phenomenologically put in by hand. This result would be helpful for 
understanding non-perturbative QCD along a parallel direction to the QCD sum rules. 


^This work was partly supported by the National Natural Science Foundation of China (NSFC) and the 
Australian Research Council 



I. Introduction 


The success of the potential model in explaining hadronic spectra and hadronic properties 
has been remarkabe. Especially, for the J/T and T families, various potential forms [1], in 
which both the Coulomb and confinement potentials are employed, can give results which are 
reasonably consistent with the data at the charm and bottom energy scales within a certain 
error. In general, it is believed that confinement comes from the non-perturbative effect of 
QCD, but unfortunately, so far, that is an unsolved problem. 

Along another line, years ago, Shifman et al. 0 proposed to study hadronic properties 
in terms of the QCD sum rules where a few non-vanishing vacuum condensates of quarks 
and gluons niq < >, ^ < > etc. describe the nonperturbative effect. They 

started from short distance, where the quark-gluon dynamics is essentially perturbative, and 
extrapolated the dynamics to larger distance by introducing non-perturbative effects step by 
step [3]. The applications of the theory extrapolated in studying hadronic properties such 
as the spectrum, decay width and hadronic matrix elements etc., indicate that one can trust 
the validity of this approach. 

Inspired by these successes, we have been trying to introduce the non-perturbative QCD 
effects, characterized by non-vanishing quark and gluon vacuum condensates, into the tradi¬ 
tional potential model [4], so that a deeper understanding of the hadronic structure and the 
underlying mechanisms which determine how quarks are bound into hadrons can be obtained. 
In our derivation of correction terms, the vacuum condensates were employed to modify the 
free gluon propagator. As a consequence, the quark-quark potential and the spectrum of 
heavy quarkonium are affected. Studying these possible effects is the main purpose of this 
paper. Meanwhile, we have noticed that beside the modification of the gluon propagator, 
characterized by the vaccum condensates, the closed-loop correction can also contribute a 
comparable effect, because these two kinds of corrections are in the same order of a*. This 
has been shown by Gupta et al [5,6], Fulcher [7] and Pantaleone et al [8]. In this investigation, 
as a preliminary study, we tempararily treat the condensate effect alone and will collect the 
closed-loop correction in our later paper. 

As to the condensate correction, there are two different results, ours [4] and Larsson’s 
[9]. The effect of the difference may show up in the spin splittings of heavy quarkonia. Thus 
our additional effort in this paper will be donated to search the possible effects of these two 
approaches. We hope this investigation may help us to understand the condensate correction 
deeply. 

As mentioned above, the non-vanishing vacuum condensates can only be used to describe 
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an extrapolation from the short distance to the medium range [3], therefore, the longer 
distance effect at the energy scale < Aqcd cannot be included in the scenario. In other 
words, this scenario is only valid within the intermediate range if only a finite number of 
vacuum condensates are kept, and for the bound states, the potential term responsible for the 
confinement should come from the larger distancel/ > Aqcd- Therefore, we are attempting 
to introduce a reasonable picture where the main contribution to the confinement is caused by 
the interaction at < Aqcd, where the physical picture is not clear yet. Since this part is not 
derivable at the present stage, we keep a phenomenological confinement form, generally the 
linear Kr term which is the main part of the confinement potential and universal to all of the 
heavy flavors. Then we perturbatively introduce the corrections induced by the non-trivial 
vacuum condensates into our framework. It would make observable contributions to some 
hadronic properties of J/T and T families, for example the spin splitting between and 
could be one of the sensitive quantities for the correction. Moreover, we would claim 
that there is no double-counting between our derived correction and the contribution from 
the larger distance 1/Aqcd- 

It should be emphasized that the amplitude of the wave function at the origin, which 
is essential to the spin splitting, depends on the potential model adopted. For instance, 
Richardson potential [10] wave function at the origin is almost the half of Cornell’s. In 
this investigation, we start with the Cornell potential, which has the simplest form among 
existant potential models, as the zeroth order approximation and then add in the condensate 
corrections in a perturbative way. Therefore, our wave function at the origin is close to that 
of the Cornell potential, and the difference between ours and Cornell’s starts to show up only 
in the approximation higher than the zeroth order. In our later work, we will study these 
differences and effects caused by adopting different model wave function as the zeroth order 
wave function. 

In the next section, the potential corrections derived with vacuum condensates are briefly 
reviewed and different formulae are analyzed. In section III, the numerical results are pre¬ 
sented and compared with the data. Finally, our results are discussed and conclusions are 
drawn. 

II. Our Model 

As SVZ suggested, there are non-zero vacuum expectation values of the quark (< 
VjgV’g >) and gluon < GG >) fields, so the propagator of the gluon should include the 
effects of these condensates. In fact, in the propagator, all the terms associated with vacuum 
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condensates are proportional to a” (n > 1). On the other hand, the contributions of higher 
order perturbative QCD corrections to the potential were discussed in Ref. [5]. In comparison 
with such corrections, the condensate terms do not suffer from the loop suppression. SVZ 
noticed that fact and then suggested that the condensates could be considered as a larger 
contribution. Therefore, at the order of we only include the condensates, but not the 
perturbative correction. 

There have been various ways to modify the potential. Except the higher order pertur¬ 
bative QCD correction, all of them attempt to include some non-perturbative effects to the 
potential, because it is sure that such effects must be taken into account. Richardson et al 
proposed a form [M] 

4 127r 1 1 

3(33-2iV^)^/„(l + |;) 

where is the momentum of the gluon exchanged between quarks and Nj represents the 
flavor number, while the linear confinement still remains unchanged as Kr. It gives an effective 
correction which indeed involves the non-perturbative effects. Besides, Fulcher gave 

V(r) = Ar — - -^---/(Ar) 

where /(Ar) has a very complicated integration form, which can be found in ref.||^]. Recently, 
Gupta et al. considered not only the higher order perturbative radiative corrections, but also 
a more complicated non-perturbative term |^]. In all these works, the corrections related to 
non-perturbative effects are phenomenologically put in according to the observation or hint 
from lattice gauge results. Some very good results which coincide with data within an error of 
a few MeV were reported 1^. We will discuss this problem in some detail in the last section. 

On the other hand, in this work we are trying to understand such corrections in terms of 
some well-established theories which can handle the non-perturbative QCD in a more natural 
way. 

Within the QCD scenario, where the non-vanishing vacuum condensates of quarks and 
gluons characterize the non-perturbative effects, the modified gluon propagator in momentum 
space can be written as [4] 
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In this expression we only keep the lowest dimensional condensates < iljq'4>q > and < GG >. 
We derived this expression in the standard way [4], in which the normal product operators 
such as tp'tpiO) and G^(0) have non-vanishing matrix elements in the physical vacuum, — i.e., 
< > and < G^ > are left as parameters to describe non-perturbative effects, and their 

values have already been determined in literature 0. 

In another way, by comparing the {2n + l)-point Green’s function (n is the number of 
external legs of V’ or A^) and the n-point Green’s function with the insertion of the operators 
V’'!/’(0) or (7^(0) (these Green’s functions are with respect to the physical vacuum), Larsson 
achieved @ 
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^ > bgl < > i(-i) _q^ 
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where the propagator was derived in the Euclidean space. Expanding the first piece with 
respect to a^, one notices that not only the coefficients at the lowest order of Og are different 
from those in Eq.(2), but also the sign of the coefficient associated with the gluon condensate 
is different from that in Eq.(2). These differences arise from different approaches used to 
deal with the non-perturbative effects and perhaps due to the the improper usage of the fixed 
point gauge which violates translational invariance. 

In this work, we also determine phenomenologically whether Eq.(2) or Eq.(3) is more 
consistent with the data. 

We write down a proper scattering amplitude between two quarks as 

M = {-igsfug{pi)-ff,Y'^q{p[)D^^’^{q'^)Uq:{p2)'yi,Y'^q'iP2) (4) 


with 

Pi - Pi = P2 -P2 = q- 

and carry out the Breit-Fermi expansion with the spinors Uq{pi) being the solutions of free 
quarks. In deriving an effective potential, the spontaneous approximation, go = 0, has been 
taken. It should be mentioned that this approximation is a traditional treatment in literatures 
although it is not quite valid in the potential derivation. Moreover, in ref.|jl^, the term g^gjy 
was kept as qiqj (go = 0). If the CVG theorem is respected, i.e. q^Uq{pi)'y^Uq{p'i) = 0, the 
term q^qu/q"^ vanishes, and then the derived potential will have a small difference with that 
without GVC theorem, even though it is not very extravagantly apart. Then one can apply 
the three-dimensional Fourier transformation to convert the propagator in momentum space 
to coordinate space and derive an effective potential between quarks. 
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The potential derived in the way used in ref. [4], in which leading order non-perturbative 
QCD effects are considered, together with the phenomenological linear confinement can be 
written in the following form: 


V{t) 


6r 


\r) + Vr^{T) + Vr'^{r), 


where is the correction from the non-trivial physical vacuum condensates, while ¥• 2 °^'' 

and are the Breit-Fermi corrections to the Coulomb and linear confinement terms, 

respectively. Due to the lengthy expressions for these potentials, we do not present them 
here. The explicit forms of these potentials can be found in Appendix. 

It is noted that: 

(1). In the preceding literature, for instance ref.[14], a comparison with the experimental 
data shows that the potential nr — ^f-, with universal values of k and Ug is applicable to 
both J/T and T in solving the Schrodinger equation. Hence, this confining potential should 
be independent of the quark mass rric or nih. In other words, it corresponds to the potential 
in the limit mg —oo (mg ^ Aqcd)- It is universal and exists in all QQ systems. On 
the other hand, all correction terms associated with vacuum condensates are of the order of 
l/rriq, which can be directly read from the Feynman diagrams. Therefore, they appear as 
mass-dependent corrections to Consequently, there is no double-counting involved in 

Kr — Moreover, the Coulomb term —comes from one-gluon-exchange and represents a 
short-distance effect where perturbative QCD works perfectly well. Therefore the correction 
induced by the condensates which manifest the non-perturbative effects does not overlap 
with the pure Coulomb part either. The nr term may be understood in the following picture. 
When mg is very large, one can define the total momentum k of the quark Q as Pq = mQV + 6 
where Q is almost on the mass shell, v is the four-velocity of Q and 6 is the so-called residual 


momentum of the order of Aqcd |13|. ^From Eq.(4), one finds that in the large mg limit the 
emitted gluons are soft. This leads to the conclusion that the confinement term, nr, plays a 
role at the energy scale Aqcd- In this picture, nr is independent of mg. 

(2) . Since all the correction terms associated with the vacuum condensates are propor¬ 
tional to 1/ rriQ , under this meaning, they are of the same order as the relativistic corrections 
in the Breit-Fermi expansion. 

(3) . Since the higher order terms are omitted, the derived potential is not appropriate in 
dealing with higher resonances. 

(4) . Our numerical results indicate that beyond the 2S state the calculated values would 
deviate from the data more and more, but for the IS, IP and 2S states, these values indeed 
make sense (see below). 
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III. Numerical Results 


In this work, to elucidate the significance of the correction, we present a few very typical 
quantities which are calculated in the framework of QCD. 


The Cornell potential —da^/Sr + nr with corresponding parameters as and k |14|, which 
gave the best fit to the J/'k and T family data, is adopted as a basic condition, and the values 
of vacuum condensates are taken from ref.|Q]. Thus there is no free parameters at all in the 
derived expressions (2) and (3). 

To our understanding, the term nr is universal to J/'k and T and dominates the con¬ 
finement part. This is consistent with the consideration in the literature which deals with 
non-perturbative corrections. Then one can consider the additional part from condensates as 
a l/rriQ correction to the potential. It is noticed that, as discussed in most of the preceding 
literatures, one always assumed that nr was caused by a scalar exchange. As a consequence, 
it would not induce a spin splitting. On the other hand, the only term responsible for the 
spin splitting is the Coulomb term which contributes as a vector potential. It should be 
mentioned that that assumption was based on phenomenological requirements — i.e. for the 
Cornell potential, if nr was induced by scalar exchange, a better fit to the spin splitting data 
could be obtained. Now, the new correction induced by the condensates contributes to the 
spin splitting, so the whole picture changes. In particular, this modification demands that 
Kr comes not only from scalar exchange, but also from vector exchange. Thus, we can write 


nr = Pur + (1 — /3)Kr 0 < /? < 1 


( 5 ) 


where the factor /3 characterizes the fraction of the confinement potential which comes from 
vector exchange, while (1 —/3) denotes that from scalar exchange. If nr is fully caused by the 
scalar exchange, then /? = 0 and it is the same assumption as in the preceding literature. The 
explicit value of /3 can be fixed by data fitting. Then the spin splitting A = 
for cc and bb systems can be calculated. The correction term which contributes to the spin 
splitting can be read as 

..2 < A “ A “>,.1 
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ref. [4] 
ref.i 


^ i < V'/?# > ref. I 
and mi = m 2 are the masses of the heavy quarks in quarkonium. 

It is noted that, the terms with in refs. [4] and Q have the same sign, but the numerical 
value of this term in ref. is three times larger than that in ref. [4]. Since the contribution from 
this term has the same sign (though < ifjxl; >) as that from the term with 6{f) (< > > 0), 

it enhances the spin splitting between 1^5o ( <a-a> = — 3 ) and ( <a-a> = 1 ). 
However, for the terms with b in refs. [4] and [P, they have not only different numerical values, 
but also opposite signs. Therefore, the term with b given by ref.j^ tends to increase the spin 
splitting, while the corresponding term in ref. [4] reduces the spin splitting. Since 

(5/288) < > r2 

r 0.4/m 


< V'/sV'/? > /mp 


0.016(—) 


Qi c 


is a small number, the existence of < GG > does not give rise to a more significant influence 
to the spin splitting than that of < >■ Anyway, a measurement of the spin splitting 

between l^S'o and tells us that the non-perturbative effects characterized by the non¬ 
vanishing vacuum condensates indeed play an important role, so manifest themselves in 
phenomenology. 

The spin splitting of l^Pj state is also influenced by the above mentioned correction as 
well as the parameter (5. The correction terms can be re-written in the following way: 



= Hr"(r) + Vf-'^(r) + + • • • 


V^corr^^) 

= H2“"(r) + V2"-°(r) + H2*^”(r) + • ■ • 

(7) 


= Vr°(r) + Vi^^(r). 
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where the superscripts cen, s — o and ten denotes the central, spin-orbit and tensor parts 
of the corresponding V)“'’^(r), respectively, and V and S represent the vector and scalar 
parts of the confinement nr in Eq.(5), respectively. Then the spectrum of l^Pj states can be 
expressed as 


M(13P2) = M + f-\g 
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M(lV) = M-f + 2g (9) 

M(l%) = M-2f-4g, 

where 

/ = < > + < V2^-°{r) > + < Vi-°{r) > 

5 = < (r) > + < yf ”(r) >, (10) 

and M is the weighted average for the l^Pj states. 

It is noticed that in Eq. dD the V and S terms are of opposite sign and the newly 
achieved correction term is opposite in sign to the Coulomb term. Then one may easily be 
convinced that to keep a good fit to data, the appearance of non-perturbative correction terms 
requires (3^0. In fact, Gupta et al. |^] also noticed that with a reasonable contribution from 
non-perturbative effects, the linear confinement nr must be a superposition of two parts as 
given in Eq.(5). As a result, they obtained an approximate /3—value of 0.25 with their model 
potential. Alternatively, with our modified potential, we have (3 ~ 0.5 ~ 0.6 for a better fit. 

The numerical calculation is performed in the following way. The Cornell potential which 
is universal to c and b families, i.e., independent of mg, is considered as the dominant part 
of the potential, and the corresponding parameters Ug and k, are determined before adding 
in the corrections. The values of ag and k are 0.381 and 0.182GeE^, respectively. And then, 
the newly derived corrections due to the non-vanishing vacuum condensates as well as the 
Breit-Fermi corrections are treated as a perturbation adding onto the dominant part. The 
resultant values for the cc and bh systems are tabulated in Table la and Ib, respectively. 
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Table la. cc system 



exp’t. 

'^Cornell _ 

■ 3 “^ + 

^(JorneLL _j_ -^corr 
_^_ycorr ycorr 

{Vr^in re/.[4]) 
d = 0.6 

■\^(JorneLL _j_ -^corr 

1 T/corr 1 T/corr 

{vr^in re/. I) 

/3 = 0.6 

'^Cornell _j_ '^corr 

'^Cornell _j_ ■^corr 
^YCorr 

P = 0.25 

l^so 

2978.8±1.9 

3074.0 

2979.3 

3026.6 

3010.8 

2999.5 

2'so 

3594.0±5.0 

3662.1 

3446.3 

3676.0 

3625.0 

3618.4 

Dsi 

3096.88±0.04 

3074.0 

3090.6 

3157.7 

3095.1 

3098.8 

2'’si 

3686.00±0.09 

3662.1 

3493.1 

3754.0 

3674.5 

3676.7 

Dpo 

3415.1±1.0 


3321.2 

3452.7 

3440.8 

3418.0 

Dpi 

3510.53±0.12 

(lPc)3497.1 

3395.1 

3533.1 

3489.3 

3480.7 

Dp2 

3556.17±0.13 


3445.1 

3605.8 

3514.4 

3527.5 

E 20 

141.07 

0 

123.9 

153.1 

73.6 

109.4 

E 21 

45.64 

0 

50.0 

72.7 

25.1 

46.7 

^ss 

118.08 

0 

111.3 

131.1 

84.3 

99.3 

^ss 

92.0 

0 

46.8 

78.0 

49.5 

58.3 


Table Ib. bb system 



exp’t. 

■^Corneii _ 

■ 3 "““ + 

■^CorneU _j_ -^corr 
_^_Y£orr ycorr 

{Vr-^in re/.[4]) 

P = 0.6 

■^O’orneff _j_ -^corr 
_^_ycorr ycorr 

{vr'^in re/.[|) 
p = 0.6 

■^Corneii _j_ -^corr 

Y^ornell _j_ -^corr 
^VCorr 

P = 0.25 

Dsi 

9460.37±0.21 

9427.0 

9466.8 

9503.8 

9452.8 

9453.7 

2'’si 

100023.30±0.31 

10007.0 

9987.7 

10073.0 

10017.0 

10017.0 

Dpo 

9859.8±1.3 


9843.6 

9907.6 

9865.8 

9861.4 

l"‘pi 

9891.9±0.7 

(lPc)9912.8 

9886.1 

9950.3 

9902.5 

9900.8 

Dp2 

9913.2±0.6 


9920.5 

9985.9 

9929.1 

9931.7 

E 20 

53.4 

0 

76.9 

78.3 

63.3 

70.3 

E 21 

21.3 

0 

34.4 

35.6 

26.6 

30.9 


In this table, E 20 = M-^3p^, E 21 = M^Sp^-M-^3p ^, and . 

The experimental data are taken from ’’Partical Physics Booklet”, July 1994, Partical Data Group. 

The decay width r(J/'I' —> e"''e“) is also a good test for various models. For r(e+e“), 
the annihilation process is related to the zero-point wavefunction of J/T. Considering the 
QCD correction, we have 


r(e'''e ) = ro(e'''e )(1 — IGas/Svr) 


and 


ro(e+e ) 


IGTregO^ 

M2 


2 


where (j){0) is the zero-point value of the J/'h wavefunction. In this scenario, the wave function 
4>{0) undergoes a modification due to the addition of the non-perturbative QCD terms. With 
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expression (2) one has 4.85X6^, while with expression (3), 4.72i^el4 and the experimental 
data is A.WKeV. Since the zero-point wavefunction is not sensitive to the new correction, the 
modified r(e^e“) does not deviate far from that calculated directly with the Cornell potential. 

IV. Discussion and Conclusion 

The potential model is successful in explaining the hadronic spectra and other properties 
of heavy quarkonia. First, people believed that quarks are confined inside hadrons, therefore 
the potential must include a confinement term. The simplest form is the linear confinement, 
kr. Besides, at short-distance, where perturbative QCD works well, the one-gluon-exchange 
provides a Coulomb-type potential. By including these two extreme sides, the Cornell po¬ 
tential, in the form of V{r) = —-|- Kr where and k are treated as free parameters, 
indeed gives reasonable results for both cc and bb families. However, there must be some 
non-perturbative effects which are not included in the linear term of the Cornell potential 
and they definitely make substantial contributions to the evaluation of spectra and other 
properties. 

In the general approaches of Richardson, Fulcher and others, the universal linear confine¬ 
ment Kr were usually kept unchanged, but the simple propagator of the gluon 
was modified by multiplying a g—dependent factor F which is model dependent. Usually, this 
factor was phenomenologically introduced based on some physical arguments or hints from 
the lattice calculations or obtained from the higher order perturbative QCD corrections Q. 

Along the other line Q, the QCD sum rules are also successful in explaining hadronic 
effects. It implies that there should be crossing between two lines. Thus the non-vanishing 
vacuum condensates which characterize the non-perturbative effects must be somewhat in¬ 
volved in the modification factor F and may be the dominant piece or at least an important 
one. 

In this work, by including the effects of quark and gluon condensates, we derive the 
modified gluon propagator. Namely, a g—dependent factor F which is similar to that shown 
in the recent literature is obtained in the framework of QCD. There are no free parameters 
in the derived expressions. 

It is important to notice that as pointed out by Shifman, this framework is an extrapolata- 
tion from short distances where perturbative QCD is reliable. Therefore, one cannot expect 
that this factor F can include as much as a purly phenomenological ansatz. But, it does shed 
light on the physical picture and enrich our understanding of the physical mechanism which 
binds quarks into hadrons. 
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For fitting experimental data, it requires that the linear confinement comes not only from 
the scalar exchange but also from the vector exchange. This is consistent with Gupta et ah’s 
model. However, the value of (3 depends on the model. In ref.[^, it is about 0.25, but in our 
case, it must be a value of 0.5~0.6, otherwise the result is not meaningful. This indicates 
that the vector exchange gives a large contribution to the linear potential. 

Our numerical results show that by considering the effects of the quark and gluon conden¬ 
sates, a more reasonable results, especially a better ht to the spin splitting E 20 , Eiq and 

( 2 ) 

can be obtained. However, does not change in the right direction. The reason is that we 
only take the lower dimensional condensates < qq > and < > into the consideration. For 

higher excited states, interaction range becomes larger, we then have to extrapolate the sce¬ 
nario to the larger distance by introducing higher dimensional condensates such as < qqG >, 

< GGG > and etc. 

There is a discrepancy between refs.[4] and |p on the sign of the coefficient of term 

< GG >. As we pointed above, maybe it is caused by using different methods or the fixed 
point gauge. Unfortunately, the contribution from this term is small compared to that from 

< qq >, therefore the numerical values calculated by using the formulae of refs. [4] and Q are 
not very far apart. We are going to pursue this problem based on both the first principle and 
phenomenology in our next work. 

In our evaluation, the closed-loop corrections are omitted since we only try to analyze 
the nonperturbative effect in the present paper. However, the closed-loop corrections are 
comparable with those from the quark and gluon condensates. Therefore, a complete analysis 
where both the condensate and loop contributions would be considered will be the aim of 
our next work. It will be helpful to clarify that the result as well as the method in refs. [4] 
is better or those in ref. [9] is more appropriate. 

It should be mentioned that in this work we take the Cornell potential as the zeroth order 
approximation. Due to the fact that the wave function at the origin depends on the model 
sensitively, it will be interesting to choose other models as the zeroth order approximation 
and compare the calculated results. This will also be done in our next work. 

Finally, the nonperturbative effects are considered by using perturbative treatment. Al¬ 
though it is simple, sometimes it does not work well. Since our purpose is to illustrate the 
effects of nonperturbative QCD effects, it would be not a serious issue. It is also interesting 
to employ other ways such as the variational method adopted by Gupta et al to re-carry out 
the analysis. It deserves further study. 
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Appendix 

By using the method mentioned in the text and ref. [4], the quark-quark potential is 
re-derived. The explicit expression of the revised quark-quark potential can be written as: 


^ ^ r(Al-Af) , . 1 , . 3 . 

U{pi,P2,f^ = -- 1 AiS{r) + A2-+ A^r + - 

4yr 4 r ^ 


+ 
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1 1 —1— /oT* I 

+ ^4 + A,r + ^ + 

1^0 


S 12 


A. 


,4 + ^ 10 - + 


fi ■ (r X pi) - /2 • (fxp 2 ) 
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where the summation over f3 is implied and the expressions of and B' are given in the 
text. 
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